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ABSTRACT 

Context. White dwarfs are nowadays routinely used as reliable cosmochronometers, allowing to date several stellar populations. 
Aims. We present new white dwarf evolutionary sequences for low-metallicity progenitors. This is motivated by the recent finding 
that residual H burning in low-mass white dwarfs resulting from Z = 0.0001 progenitors is the main energy source over a significant 
part of their evolution. 

Methods. White dwarf sequences have been derived from full evolutionary calculations that take into account the entire history of 
progenitor stars, including the thermally-pulsing and the post-asymptotic giant branch phases. 

Results. We show that for progenitor metallicities in the range 0.00003 < Z < 0.001, and in the absence of carbon enrichment due 
to the occurrence of a third dredge-up episode, the resulting H envelope of the low-mass white dwarfs is thick enough to make stable 
H burning the most important energy source even at low luminosities. This has a significant impact on white dwarf cooling times. 
This result is independent of the adopted mass-loss rate during the thermally-pulsing and post-AGB phases, and the planetary nebulae 
stage. 

Conclusions. We conclude that in the absence of third dredge-up episodes, a significant part of the evolution of low-mass white dwarfs 
resulting from low-metallicity progenitors is dominated by stable H burning. Our study opens the possibility of using the observed 
white dwarf luminosity function of low-metallicity globular clusters to constrain the efficiency of third dredge up episodes during the 
thermally-pulsing AGB phase of low-metallicity progenitors. 
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1. Introduction 


White dwarf stars are the most common end-point of stellar evo¬ 
lution. These old stellar remnants encode valuable information 
about the evolutionary history of their progenitors, providing 
a wealth of information about the physical evolutionary pro¬ 
cesses of stars, the star formation history, and about the char¬ 
acteristics of various stellar populations. Moreover, their struc¬ 
ture and_eyolutionary_propertes are now well understoo d 


Fontaine & B rassardI (l2008l) : IWinget & Kenled (l2008h . and 
lAlthaus et al. ( 2010bl) for specific reviews. Hence, white dwarf 


cooling ages are currently considered one of the best age indi- 


and globular clusters fseelWinget et al.l2009 

^ ^ . 

: Garcfa-Berro et al.l 

120 lOl 

Jefferv et alj|201 li Bono et aLll2013l: 

Hansen et al. 2013, 


for some recent applications). Consequently, the use of white 
dwarfs as reliable clocks to date stellar populations has led 
to renewed efforts in computing full evolutionary models for 
these stars, taking into account all the relevant sources and 
sinks of energy, and the evolutionary h istory of progenitor stars 
(iRenedo et al.luOTotlSalaris et aLll^lOl) . 


It has been generally accepted that white dwarf evolution is 
almost insensitive to the initial metal content of progenitor star. 
This is a consequence of the strong gravity that characterizes 


white dwarf atmospheres. The effect of a large surface grav¬ 
ity is that metals in the very outer layers of these stars are de¬ 
pleted in relatively short timescales, and thus the cooling rates 
do not depend on t he initial metal content - see, for instance, 
ISalaris et al.l d2010ll . Accordingly, age determinations based on 
white dwarf cooling times are not largely affected by uncertain¬ 
ties in the determination of the metallicity of the parent popula¬ 
tion. However, the metallicity of the progenitor stars may play 
some role in determining the core chemical composition of a 
given white dwarf, and thus could affect the cooling rates. In 
particular, because of difference in the initial progenitor mass - 
that could result from the dependence of the initial-to-final-mass 
relation on the progenitor metallicity - different oxygen profiles 
for a given white dwarf mass could be expected depending on 
the progenitor metallicity. 

Nevertheless, the impact of progenitor metallicity on white 
dwarf evolution c ould be much more re l evant than hitherto as¬ 
sumed. Recently, iMiller Bertolami et al.l (l2013ll have computed 
for the first time white dwarf sequences resulting from the full 
evolution of progenitor stars with very low metallicity (Z = 
0.0001). They found that due to the low metallicity of the pro¬ 
genitor stars, white dwarfs are born with thicker H envelopes, 
leading to more intense H burning shells as compared with 
their solar metallicity counterparts. Specifically, they found that 
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Table 1. Basic model properties for sequences with Z = 0.00003, 0.0001, 0.0005 and 0.001. 


A/zams (A/q) 

A/wd (Mq) 

log Mh (Mq) tus (Gyr) 

tHe (Gyr) 

tiTp (Gyr) 

A^tp 

C/O 

Z = 0.00003 

0.80 

0.50611 

-4.253 

11.828 

12.602 

12.701 

3 

16.221 

0.90 

0.52007 

-3.983 

7.741 

8.312 

8.403 

3 

10.618 

0.95 

0.53487 

-3.859 

6.388 

6.888 

6.978 

2 

31.084 

1.00 

0.53967 

-3.982 

5.330 

5.771 

5.859 

3 

20.180 

1.50 

0.61840 

-3.949 

1.304 

1.485 

1.623 

12 

23.882 

1.75 

0.68029 

-4.224 

0.773 

0.903 

1.042 

19 

23.647 

2.00 

0.77021 

-4.376 

0.504 

0.586 

0.711 

25 

12.736 

Z = 0.0001 

0.80 

0.51976 

-3.220 

11.863 

12.741 

12.844 

1 

0.270 

0.85 

0.53512 

-3.304 

9.500 

10.269 

10.367 

3 

0.301 

0.90 

0.54839 

-3.360 

7.694 

8.347 

8.440 

3 

0.295 

0.95 

0.56145 

-3.435 

6.329 

6.901 

6.998 

4 

0.324 

1.00 

0.56765 

-3.461 

5.269 

5.793 

5.886 

4 

0.332 

1.25 

0.61940 

-3.652 

2.462 

2.793 

2.886 

8 

0.325 

1.50 

0.66588 

-3.851 

1.276 

1.486 

1.580 

12 

0.264 

2.00 

0.73821 

-4.347 

0.581 

0.605 

0.749 

23 

11.477 

2.50 

0.82623 

-4.648 

0.342 

0.348 

0.421 

35 

13.355 

Z = 0.0005 

0.80 

0.50932 

-3.231 

12.008 

13.124 

13.242 

1 

0.241 

0.85 

0.54164 

-3.424 

9.551 

10.520 

10.622 

3 

0.296 

0.90 

0.54619 

-3.432 

7.801 

8.601 

8.706 

3 

0.240 

1.00 

0.56634 

-3.512 

5.335 

5.955 

6.059 

5 

0.267 

1.25 

0.60391 

-3.717 

2.392 

2.797 

2.898 

7 

0.304 

1.50 

0.65406 

-3.904 

1.337 

1.539 

1.638 

12 

0.298 

2.00 

0.71244 

-4.245 

0.610 

0.633 

0.799 

20 

3.202 

2.50 

0.81197 

-4.485 

0.357 

0.364 

0.441 

31 

1.289 

Z = 0.001 

0.85 

0.53846 

-3.434 

9.892 

11.033 

11.136 

2 

0.307 

1.00 

0.55946 

-3.561 

5.411 

6.210 

6.313 

4 

0.299 

1.25 

0.60195 

-3.747 

2.384 

2.906 

3.007 

7 

0.289 

1.50 

0.63962 

-3.924 

1.386 

1.590 

1.687 

11 

0.278 

1.75 

0.66940 

-4.035 

0.898 

0.983 

1.092 

15 

0.269 

2.25 

0.75394 

-4.329 

0.477 

0.489 

0.605 

26 

0.262 


Notes. Mzams: initial mass, Mwd: white dwarf mass, log Mh: logarithm of the mass of H at the maximum effective temperature at the beginning 
of the cooling branch, Jms^ age at the end of the main sequence (defined at the moment when the central H abundance becomes 10“®), age at 
the beginning of the core He burning, fixp: age at the first termal pulse, Njp: number of thermal pulses, C/O: surface carbon to oxygen ratio after 
departure from the AGB, at log Teff - 4. 


for this metallicity and for white dwarf masses smaller than 
~ 0.6 Mq, nuclear reactions are the main contributor to the stel¬ 
lar luminosity for luminosities as low as log(L/LQ) ^ -3.2, de¬ 
laying the white dwarf cooling for significant time intervals. This 
is in contrast with the standard assumption of neglecting residual 
H burning made in all the existing calculations of white dwarf 
cooling. This assumption is usually well justified because sta¬ 
ble H shell burning is expected to be a minor source of energy 
for stellar luminosities below ~ 100 Lq. However, H burning 
never ceases completely, and depending on the mass of the white 
dwarf and on the precise mass of H left during the previous evo- 
lutionary phases - which depends critically on metallicity, see 
llben & MacDonald! (1 19861) - it may become a non-negligible en- 
ergy source for white dwarfs with H atmospheres. As shown by 
iMiller Bertolami et al.l (l2013h a correct assessment of the role 
played by residual H burning during the cooling phase requires 
a detailed calculation of the white dwarf progenitor history, par¬ 
ticularly during the thermally-pulsing asymptotic giant branch 
(AGB) phase. 

In view of this important finding, here we present new white 
dwarf evolutionary sequences for a wide range of metallicities 
that cover the metal content expected in some metal-poor stel¬ 
lar populations, like the galactic halo or the Galactic popula¬ 


tion of globular clusters. These sequences are of interest for the 
cosmochronology of old stellar systems. We compute the full 
evolution of white dwarf stars by consistently evolving the pro¬ 
genitor stars through all the stellar phases, including the ZAMS, 
the red giant branch, the core helium flash (whenever it occurs), 
the stable core helium burning, the AGB phase and the entire 
thermally-pulsing and post-AGB phases. With the aim of deter¬ 
mining at which metallicity residual H burning becomes relevant 
for white dwarf evolution, we have considered four progenitor 
metallicities: Z = 0.001, 0.0005, 0.0001, and 0.00003. To the 
best of our knowledge, this is the first set of fully evolution¬ 
ary calculations of white dwarfs resulting from low-metallicity 
progenitors covering the relevant range of initial main sequence 
and, correspondingly, white dwarf masses. We emphasize that 
our complete evolutionary calculations of the history of the pro¬ 
genitors of white dwarfs allowed us to have self-consistent white 
dwarf initial models. That means that in our sequences the resid¬ 
ual masses of the H-rich envelopes and of the helium shells were 
obtained from evolutionary calculations, instead of using typical 
values and artificial initial white dwarf models. 

The paper is organized as follows. In Sect. |2]we briefly de¬ 
scribe our numerical tools and the main ingredients of the evolu¬ 
tionary sequences, while in Sect. [3] we present in detail our evo- 


2 



















Althaus et al.: White dwarf evolutionary sequences for low-metallicity progenitors 


lutionary results. There, we also explore the dependence of the 
hnal H content on the mass loss and the treatment of overshoot¬ 
ing during the thermally-pulsing AGB phase. Finally, in Sect. |4] 
we summarize the main hndings of the paper, and we elaborate 
our conclusions. 


2. Numerical setup and input physics 

The evolutionary calculations presented in this w ork have been 
done using the LPC ODE stellar evolutionary code (lAlthaus et alJ 
I 2 OO 3 LI 2 OO 5 L l2012h . This is a well-tested and callibrated code 
that has been amply used in the study of different aspects 
of low-mass star evolution. In particular, it has been em¬ 
ployed to compute very accurate and realistic white dwarf 
models, including the formati on and evolution of ex t remel y 

iMiller Bertolami et aTI (l2008l): 


low-mass white dwarfs - 
iGarcfa-B erro et aTI (l2010t); Althaus eTakri loiOal): Rened^o et al.l 


(I2OIOI): iMiller Bertolami et al.l ( ^1 Ih: IWachlin et alJ ( 2011 ): 

ICorsico et al.l ~( 20121) and lAlthaus et alJ ( 2013h and references 
therein for recent applications. LPCODE has been tested against 
another white dwarf evolutionary code, and uncertainties in the 
white dwarf cooling ages arising from the different numerical 
implementations o f the stellar evoluti on equations were found 
to be below 2% (ISalaris et al.l l2013h . A detailed description 
of the input physics and numerical procedures can be found 
in these works. Briefly, extra mixing due to diffusive convec¬ 
tive overshooting has been considered during the core H and 
He burning, but not during the thermally-pulsing AGB. The 
breathing pulse instability occurring towards the end of core 
helium burning has been suppressed. The nuclear network ac¬ 
counts explicitly for the following elements: *H, ^H, ^He, "^He, 
^Li, ^Be, '2c, '3C, '“'N, '^N, '"O, '^O, '*0, '“'F, ^ONe and 
22Ne, together with 34 thermonuclear reaction rates for the pp- 
chains, CNO bi-cy cle, and helium burn ing that are identical to 
those described in lAlthaus et al.l (l2005l) . with the exception of 
the reactions '^C + p^ '^N - H y —> '^C + e^ + Vg and 
*2C(p,y)*"'N, which are taken from lAngulo et al.l(ll999ll. In addi¬ 
tion, th e reacion rate *"'N(p,y)*^0 was taken from llmbriani et aP 
( 20051). Radiative and con ductive o pacities are taken fro m OPAL 
( Iglesias & Rogerslll996h and from ICassisi et al.l (l2007h . respec¬ 
tively. The equation of state during the main sequence evolu¬ 
tion is that of OPAL for H- and helium-rich composition, and a 
given metallicity. Updated low-temperature molecular opacities 
with varying C/O ratios are used, which is relevant for a realistic 
treatment o f progenitor evolution du ring the thermally-pulsing 
AGB phase dWeiss & Fergusonl2009l) . Fo r this purpose, we have 
adopted the low- temperature opacities of Ferguson et al.l (l2005l) 
and presented in IWeiss & Fergusoril (l2009l) . In LPCODE molecu- 
lar opacities are computed by adopting the opacity tables with 
the correct abundances of the unenhanced metals (e.g., Fe) and 
C/O ratio. Interpolation is carried out by means of separate 
cuadratic interpolations in /? = p/Tg^, T and Ah, but linearly 

in Ac/Ao- 

For the white dwarf regi me, we use the equation of state 
of iMagni & Mazzitellil (Il979l) for the low-density regime, while 
for the high-density re gime we consider the equation of state of 
ISegretain et al.l (119941) . which accounts for all the important con¬ 
tributions for both the liquid and solid phases. We also take into 
accou nt the effects of element diffusion due to gr avitational set¬ 
tling (llsern et al.ll 19911: iGarcfa-Berro et al.ll2008l) . chemical and 
thermal diffusion of *H, ^He, "'He, '^C, '^C, '"'N and '®0, see 
lAlthaus et al.l (l2003l) for details. For the white dwarf regime and 
for effective temperatures lower than 10,000 K, outer bound¬ 
ary conditions for the evolving models are derived from non¬ 


grey model atmospheres (iRohrmann et al.l [20121) . All relevant 
energy sources are taken into account in the simulations, in¬ 
cluding nuclear burning, the release of latent heat of crystal¬ 
lization, as well as the release of the gravitational energy as¬ 
sociated with the carbon-oxygen phase separation induced by 
crystallization. In particular, H burning is considered down to 
log(L/L 0 ) -4. The inclusion of all these energy sources is 
done self-consistently and loc ally coupled to the fu ll set of equa¬ 
tions of stellar evolution - see lAlthaus et al.l (1201 2l) for a detailed 
description. 

The full calculation of the evolutionary stages leading to the 
formation of the white dwarf is mandatory for a correct eval¬ 
uation of the possible role of residual nuclear burning in cool 
white dwarfs. Th e present vers i on of L PCODE improves consider¬ 
ably that used in iRenedo et al.l (l2010l) . both in the microphysics, 
and the numerical treatment of stellar equations and mesh dis¬ 
tribution. This has allowed us to follow the progenitor evolution 
through the entire thermally-pulsing AGB phase. In fact, in con¬ 
trast with this previous study, in which was necessary to artifi¬ 
cially modify the opacity profile close to the base of the convec¬ 
tive envelope to circumvent the numerical difficulties arising at 
tfie end of the thermally-pulsing AGB, in the present work all 
of the progenitors have been evolved naturally throughout the 
entire thermally-pulsing AGB, that is, we have not forced our 
models to abandon early the thermally-pulsing AGB phase in 
order to get realistic white dwarf structures. 

In computing our new set of white dwarf sequences, over¬ 
shooting was disregarded during the thermally-pulsing AGB 
phase. Our reason for this choice is that we want to constrain 
how important can residual H burning be for the evolution of 
low-mass low-metallicity white dwarfs. As a result, it is ex¬ 
pected that the mass of the H-free core of our sequences grad¬ 
ually grows as evolution goes through the TP-AGB. In fact, 
the inclusion of overshooting at all convective boundaries dur¬ 
ing the TP-AGB leads to a strong enhancement of third dredge 
up events, which p revent the growth of the hydrogen-free core 
(ISalaris et al.l2009h . In the absence of overshooting, third dredge 
up in low mass stars in very feeble. Hence, mass loss plays a 
major role in determining the final mass of the H-free core at 
the end of the TP-AGB evolution, and thus the initial-final mass 
relation. In tfiis work, ma ss loss during tfie RGB was taken from 
|Schroder_&Qmtz (l200-5h. Fo r tfie AGB pfiase, we use again tfiat 
of Scfiroder & Cuntzl (l2005l) for pulsation periods shorter than 


50 days. For l onger periods, mas s loss i s tak en as the maximum 
of the rates of I Schroder & Cuntzl (l2005|) and iGroenewegen et al.l 
(I 2 OO 9 I) for oxygen r i ch sta rs, or the maximum of the r ates 
of ISchroder & (ihint3 (l2005l) and iGroenewegen et aP (Il998l) for 
carbon-rich stars. In all of our calculations, mass loss was sup¬ 
pressed after the post-AGB remnants reach log T^ff - 4. This 
choice does not bear consequences for the final mass of the H 
content that will characterize the white dwarf, see next section. 

The initial models for our white dwarf sequences correspond 
to stellar configurations derived from tfie full evolutionary calcu¬ 
lations of tfieir progenitor stars. Four metallicities for tfie progen¬ 
itor stars fiave been considered: Z = 0.001,0.0005,0,0001 and 
0.00003. The initial H abundances, by mass, are, respectively 
0.752,0.7535,0.7547, and 0.7549. A total of 30 evolutionary se¬ 
quences have been computed self-consistently from the ZAMS, 
through the giant phase and helium core flash, the thermally- 
pulsing phase, and the planetary nebulae stage. They encompass 
a range of initial stellar masses from 0.80 to 2.5 Mq. The initial 
mass of the progenitor stars at the ZAMS and the resulting white 
dwarf mass, together with other quantities that will be discussed 
below, are fisted in Table [1] In all cases, the white dwarf evolu- 


3 











































































































Althaus et al.: White dwarf evolutionary sequences for low-metallicity progenitors 



^wd^^sun 

Fig. 1. Logarithm of the total mass of H (in solar masses) at the 
point of maximum effective temperature of the cooling sequence 
versus the white dwarf mass for Z = 0.0001,0.0005 and 0.001. 
In addition, we show the predictions of iRenedo et alJ (l2010h for 
Z = 0.001 and Z = 0.01 progenitors. 


tion has been computed down to log(L/LQ) = -5.0. Moreover, 
additional sequences have been computed to assess the impact 
on the final H content due to mass loss and overshooting (and 
the occurrence of third dredge-up in low-mass progenitors) dur¬ 
ing the thermally-pulsing AGB phase, see Sect. 13.11 

3. Evolutionary results 

We will see that in the absence of third dredge-up episodes dur¬ 
ing the AGB phase, most of the evolution of white dwarfs result¬ 
ing from low-mass (M ^ 1.25 Mq) low-metallicity progenitors 
is dominated by stable H burning. In particular, the occurrence of 
third dredge-up reduces the hnal mass H content with which the 
white dwarfs enter their cooling track, thus weakening the role 
played by residual H burning (see Sect. [IB- Contrary to the 
situation of metal-rich progenitors, for which there exist direct 
observational inferences about the occurrence of third dredge- 
up in AGB stars with stellar masses larger that M ~ 5 Mq 

dLebzelter et alJl^OSl) . the situation for their metal-poor coun¬ 
terparts seems to be not clear, particularly for the low-mass AGB 
stars. 

Evidence pointing towards the occurrence of third dredge-up 
in progenitor stars with very low metallicity is provided by ob¬ 
servations of carbon-enhanced metal-poor stars (CEMP) stars. 
However, the origin of the carbon-enhancement for CEMP stars 
is not fully understood. A possible scenario involves binary pair¬ 
ings of a low-mass star (the one observed today) and a more 
massive binary companion that can donate carbon in a mass- 
transfer event. This requires that the more massive star to be in 
the right mass range to pass through an AGB phase with sig¬ 
nificant production of car bon (and s-process elements), roughly 
between 1.3 and 6.0 Mq (ICarollo et al.ll2014ll . Thus, the CEMP 
phenomenon is associated with massive and intermediate-mass 
progenitor stars, and not with the low-mass stars that are relevant 
for our study. It is worth noting that the third dredge-up evidence 



Fig. 2. Ratio of the carbon to oxygen abundances at the end of the 
thermally-pulsing AGB phase as a function of the white dwarf 
mass for Z = 0.0001,0.0005, and Z = 0.001 progenitors - solid, 
dotted and dashed lines, respectively. In addition, we show with 
hlled and hollow squares the predictions from our low-mass Z = 
0.0001 progenitors in the case that overshooting is considered 
during th e thermally-pulsing AGB phase and the corresponding 
results of lWeiss & FergusonI (l2009l) for Z = 0.0005 progenitors. 


suggested by the CEMP stars is in line with our calculations for 
the Z = 0.0001 progenitors, which predict a strong carbon en¬ 
richment even in the absence of overshooting for stellar masses 
larger than ~ 1.5 Mq (Table[B. 

This evidence is not conclusive about the occurrence of third 
dredge-up episodes in low-mass, low-metallicity AGB stars. On 
the theoretical side, turning on convective overshooting at the 
base of the convective en velope artihcial ly triggers the occur¬ 
rence of third dredge-up (lMowlavilll999h . But the initial-final 
mass relation of solar metallicity stars suggests that third dredge- 
up episodes should be less relevant than predicted by evolution¬ 
ary models that include oversho oting at all convecti ve bound¬ 
aries during the TP-AGB phase dSalaris et al.ll200^ . This ex¬ 
pectation is in line with simple hydrodynamical arguments sug¬ 
gesting that the penetration of the convective envelope though 
the H/He transition is expected to be minor. A simple estimate 
can be done on the basis of the bulk R ichardson number (Ris). 
In particular, iMeakin & Arnett! (l2007l) proposed that the bulk 
Richardson number could be used as a dimensionless estimation 
of the “stiffness" of a stable layer against convection. Thus, this 
dimensionless number could be used to constrain the strength of 
extra-mixing episodes. These authors suggest that high values 
of RIb should inhibit the occurrence of apreciable extra-mixing. 
This is in line with terrestrial simulations which show that the 
entrainment proc ess in the ocean is estim ated to operate up to 
Rie ~ 10^ - 10® (IMeakin & Arnetlll200^ . 

In order to assess the possibility that the outer convection 
zone penetrates into the carbon-rich layers, we have calculated 
RIb along de H/He transition f or our 0.90Mn sequence with 
Z = 0.0001, following Eq. (2) in IMeakin & Arnett (120071) . This 
calculation was performed at the moment when the outer con- 


4 






































Althaus et al.: White dwarf evolutionary sequences for low-metallicity progenitors 



Log(LXun) 



Log (L^/Lguj^) 



LogCL'X 



Q ..I.. I I i d 

1 0 - 1-2 -3 -4 -5 


Log {L*/Lgup) 


Fig. 3. Fraction of the total white dwarf luminosity due to H nuclear reactions for the four different metallicities studied in this work. 
Note that for Z = 0.0001 and Z = 0.0005, residual H burning becomes the main energy source in low-mass white dwarfs, even at 
low luminosities. 


vective zone reaches its maximum depth, during the second ther¬ 
mal pulse, when third dredge-up episodes are expected to occur. 
We hnd that the Rie values range from ~ 1 in the H rich zone, 
up to ~ 5 X 10^ at the outermost layers that has been enriched in 
carbon by the pulse driven convection zone. Although this does 
not mean that some extra-mixing should not occur, this last value 
is high enough to prompt us to suspect that is not entirely cor¬ 
rect to assume the occurrence of extra-mixing through the entire 
H-He transition layers. Hence, on the basis of this hydrodynam- 
ical argument, it cannot be conclusive about the occurrence of 
appreciable third dredge-up episodes in these low-mass stars. In 
view of these results, we believe to be worthwhile to compute 
sequences considering and disregarding extra-mixing during the 
thermally-pulsing AGB phase, and explore the consequences for 
white dwarf evolution. 

In all of our sequences, evolution starts from the ZAMS and 
is followed through the stages of stable H and helium core burn¬ 
ing, the stage of mass loss during the entire thermally-pulsing 
AGB, the domain of the planetary nebulae at high effective tem¬ 
perature, and hnally the terminal white dwarf cooling track, un¬ 
til very low surface luminosities are reached. For most of them, 
evolution has proceeded to the helium core flash on the tip of the 
Red Giant Branch, and the following recurrent sub-flashes, be¬ 
fore the progenitors reach the stable core helium-burning stage 
on the Horizontal Branch. The final white dwarf mass (in solar 


masses) of these sequences can be found in Table [T] together 
with the corresponding evolutionary times (in Gyr) for the main 
relevant stages of progenitor stars. Also listed in this table are 
the total mass of residual H at the point of maximum effective 
temperature at the onset of the cooling track, the number of ther¬ 
mal pulses on the AGB, and the final C/O ratio after the end of 
the thermally-pulsing AGB phase. 

An expected result shown in Table [T] is that the residual H 
content left by prior evolution decreases with increasing white 
dwarf masses, a fact that helps to understand the dependence 
of residual nuclear burning on the stellar mass discussed below. 
Another trend that can be seen in Table [1] is that the residual 
H content displays a marked dependence on the metallicity of 
the progenitor. This is related to the well-known result that low- 
metallicity progenitors depart from th e AGB with more massive 
envelopes dlben & MacDonald 1 1986h . leading to white dwarfs 
with thicker H envelopes. Indeed, we And that, in general, for a 
given white dwarf mass, a larger H content is expected to be left 
as the metallicity progenitor decreases, the only exception being 
the sequences with the lowest metallicity (Z = 0.00003). This 
can best appreciated in Fig.[Tl where the total mass of H (in solar 
units) obtained in our calculations at the beginning of the cooling 
branch is shown for Z = 0.0001,0.0005 and 0.001 as a function 
of the white dwarf stellar mass. In addition, the predictions of 
iRenedo et al.l (l2010l) for Z = 0.001 and Z = 0.01 progenitors are 
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shown. As a result of the more massive H envelopes, residual H 
burning is expected to become more relevant in low-mass white 
dwarfs with low-metallicity progenitors, see later in this section. 

We will show in Sect. 13.11 that for a given progenitor 
metallicity, the amount of H with which the remnant will en¬ 
ter the cooling branch is intimately connected with the oc¬ 
currence of the third dredge-up during the thermally-pulsing 
AGB phase. During the third dredge-up, the envelope becomes 
carbon-enriched as a result of surface convection penetrating 
into carbon-rich deeper layers. In this regard in Fig. |2]we dis¬ 
play the final ratio of carbon to oxygen abundances at the end 
of the thermally-pulsing AGB phase, as a function of the white 
dwarf mass for Z = 0.0001,0.0005, and 0.001. This figure shows 
that for lower metallicities the effect of the third dredge-up on 
the C/O ratio is pronounced. We note that in our calculations 
we did not find any third dredge-up episode except for the two 
more massive sequences - those with initial ZAMS masses 2.0 
and 2.5 Mq - with Z = 0.0001 and Z = 0.0005. The occur¬ 
rence of third dredge-up in these sequences is clear from Figl2] 
where a sudden increase in the hnal C/O ratio can be seen - at a 
white dwarf mass of ~ O.65M0. The third dredge-up reduces the 
amount of mass of residual H, as it is inferred from the change in 
the corresponding curves shown in Fig[T] Since the role played 
by nuclear burning during white dwarf evolution is strongly tied 
to the residual H mass, which is essentially a function of the 
occurrence of third dredge-up and the metallicity of the progeni¬ 
tor, our derived values for the H envelopes should be considered 
as upper limits. In fact, as already mentioned, we disregarded 
overshooting during the thermally-pulsing AGB phase. In par¬ 
ticular, in our simulations third dredge-up is much less efficient 
in the absence of overshooting. Thus our choice of neglecting 
overshooting favors the formation of white dwarfs with larger 
amounts of H (see Sect l3.1l) . 

In connection with this, we mention that the mass of the H 
that is left on the white dwarf (listed in Table [T]) is indepen¬ 
dent on the occurrence of mass loss during the planetary nebulae 
stage. As mentioned before, for all our sequences mass loss was 
suppressed after the post-AGB remnants reach log - 4. To 
test the impact of additional mass loss on the residual mass of H, 
we have re-computed the post-AGB evolution of our 0.80 Mq 
sequence with Z = 0.0005 from this effective temperature, 
but now considering mass loss during all the planetary nebulae 
stage, including the hot stages of white dwarf evolution, follow¬ 
ing the prescription of the mass loss rates for hot central stars 
of plan etary nebulae provided by iBloeckeil (Il995h from the re¬ 
sults of iPauldrach et al.l (Il988h . We find that the mass of H with 
which the white dwarf enters its cooling track is not altered by 
the occurrence of mass loss during these stages. Mass loss dur¬ 
ing this phase increases the rate of evolution of the post-AGB 
remnant (since the H envelope is reduced now by both nuclear 
burning and mass loss) with the result that the remnant reaches a 
given effective temperature with the same mass of the H content 
irrespective of whether mass loss is considered or not. 

We remark that the residual masses of H in our Z = 0.001 
seque nces are somewhat la rger than those found in our previous 
study (iRenedo et al.ll201^ for the same metallicity, see Fig. [1] 
This might be due to the several numerical and physical updates 
in LPCODE. In particular, in contrast with our previous study, here 
we have not artificially altered the opacity at the base of the con¬ 
vective envelope in order to get initial white dwarf structures (see 
Sect.|2]l. Instead, our sequences were computed self-consistently 
throughout this stage and during the departure from the AGB. 

As a result of the more massive H envelopes, residual 
H burning is expected to become more relevant in white 



^coolingl^y^] [Myr] 


Fig. 4. Impact of residual nuclear burning on the cooling times of 
three selected white dwarf sequences resulting from progenitors 
with Z = 0.0001. The solid line illustrates the evolution in the 
case in which FI burning is taken into account, whereas the dot¬ 
ted line shows the case in which nuclear burning is disregarded. 
For the 0.53512 Mq white dwarf cooling sequence we also show 
the evolution when diffusion is neglected (dashed line). 


dwarfs with low-metal l icity p rogenitors. This was borne out by 
iMiller Bertolami et alJ d2013l) . who found that for white dwarfs 
with progenitors of metallicity Z = 0.0001 stable H burning 
becomes one of the main energy sources of low-mass white 
dwarfs for substantial periods of time. The calculations per¬ 
formed here show that this is also true even for white dwarfs re¬ 
sulting from progenitors with metallicities as large as Z » 0.001. 
This is inferred from Fig. |2 where we plot the fraction of the 
total luminosity that results from H burning at different stages 
of the white dwarf cooling phase for the four metallicities ex¬ 
plored in our work. It is clear that the luminosity of low-mass 
white dwarfs resulting from progenitors with metallicity in the 
range 0.0001 $ Z ^ 0.001 is completely dominated by nuclear 
burning, even at rather low luminosities. Specifically, note that 
for white dwarfs with M ^ 0.6 Mq the energy release of H 
burning is the main energy source at intermediate luminosities 
(-3.2 ^ log(L/LQ) ^ -1, or 7800 ^ T^siK) ^ 26,000). In 
sharp contrast, in the case of white dwarfs resulting from pro¬ 
genitors with Z = 0.00003, the role played by residual H burning 
is much less relevant. At such very low metallicities, all the pro¬ 
genitor sequences we computed experience a H ingestion flash 
during the first thermal pulses, resulting in a strong carbon en¬ 
richment of the outer layers (see Table [T]i, and in white dwarfs 
with thin H envelopes that are unable to sustain appreciable nu¬ 
clear burning. Hence, according to our calculations, progenitors 
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Fig. 5. Same as Fig. |4|but for the white dwarf sequences of pro¬ 
genitors with Z = 0.0005. 


with metallicities lower than Z ~ 0.00003 are not expected to 
form white dwarfs in which residual H burning is appreciable. 

The impact of residual H burning on the cooling times is 
depicted in Figs. |4| |5] and | 6 ] for progenitor metallicities of 
Z = 0.0001, 0.0005 and 0.001, respectively, and for three se¬ 
lected white dwarf sequences. The solid lines display the cooling 
times (in Gyr) of the sequence in which H burning is considered, 
while the dotted lines correspond to the situation in which resid¬ 
ual H burning is disregarded. The right panels of each figure 
display the corresponding delays (in Myr) introduced by nuclear 
burning. Clearly, nuclear burning leads to substantial delays in 
the cooling times, as compared with the situation encountered in 
white dwarfs that result from Solar metallicity progenitors, for 
which nuclear burning does not play a leading role, and most 
of the energy release comes from the thermal energy stored in 
the interior. This is most noticeable for those white dwarfs that 
resu lt from progenitor stars with metallicity of Z = 0.0001 - 
(see lMiller Bertolami et alj|2013h - but also for white dwarfs re¬ 
sulting from progenitors with larger metallicities (Z = 0.0005). 
Neglecting the energy released by nuclear burning leads to a 
marked underestimation of the cooling times, by more than a fac¬ 
tor of 2 at intermediate luminosities. Even at log(L/L 0 )-4, 

residual H burning leads to an increase in the cooling times of 
the low-mass white dwarfs by 20 - 40% (see Table|2]). Note also 
that for more massive white dwarfs, the impact of H burning 
on cooling times is much less relevant. In the case of progen¬ 
itors with larger metallicity (Z = 0.001), the delay introduced 
by H burning is also apparent in the least massive white dwarfs, 
and, as expected from the earlier di scussion, more relev ant than 
found in our previous calculations (iRenedo et alJl201^ for the 
same metallicity. Here, H burning increases the cooling times 


Fig. 6. Same as Fig.|4|but for the case of white dwarf progenitors 
with Z = 0.001. 


by 10 - 20% at log(L/LQ) -4 in the case of low-mass white 

dwarfs. It is clear that stable nuclear burning in low-mass, low- 
metallicity white dwarfs can constitute a main energy source, 
delaying substantially their cooling times at low luminosities. 

To assess the influence of element diffusion on H burning, we 
have recomputed the evolution of the 0.53512 Mq white dwarf 
sequences with Z = 0.0001 for the extreme situation of neglect¬ 
ing all of the diffusion processes. The resulting cooling time is 
shown in Fig. |4|as a dashed line. Note that even in this extreme 
case, substantial delays in the cooling times are expected as a 
result of H burning. 

Figure |7] shows the evolution of the mass of H for three se¬ 
lected white dwarf sequences with Z = 0.0001. The mass of 
the initial H layer decreases - by almost a factor of two - for 
a period during which H burning supplies most of the surface 
luminosity of the white dwarf. However, as soon as the mass of 
the H layer decreases below a certain threshold, which depends 
on the white dwarf mass, the pressure at the bottom of the en¬ 
velope is unable to no longer support nuclear reactions, and the 
total H mass of the model reaches a constant value. For com¬ 
parison purposes, the run of the mass of H for the 0.53512 Mq 
white dwarf sequence when element difussion processes are ne¬ 
glected is also shown in the figure (dashed line). In the absence 
of element diffusion, a larger fraction of the H content of the star 
burnt at higher luminosities, as compared with the case in which 
diffusion is allowed to operate. 

In Fig.jSjwe show the time dependence of the various energy 
sources and sinks during the white dwarf cooling phase for the 
0.53512 M 0 and the 0.82623 M 0 sequences. These sequences 
result, respectively, from progenitor stars with masses 0.85 and 
2.5 Mq, and a metallicity Z = 0.0001. For the more massive 
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Table 2. White dwarf cooling ages for H burning sequences and their corresponding time delays compared to the sequences in 
which nuclear burning is neglected. 


-log(L/LQ) 





tcooi (Gyr) 







Z = 0.0001 



Z = 0.0005 



Z = 0.001 



1 0.535 M© 

0.561 Mq 

0.666 Mq 

0.542 Mq 

0.566 Mq 

0.712 M© 

0.538 Mq 

0.559 Mq 

0.669 Mq 

1.0 

0.04 

0.03 

0.02 

0.03 

0.03 

0.01 

0.03 

0.02 

0.01 

2.0 

0.40 

0.33 

0.21 

0.29 

0.27 

0.19 

0.26 

0.22 

0.19 

3.0 

1.95 

1.67 

1.01 

1.54 

1.23 

0.92 

1.20 

1.07 

0.91 

3.5 

3.13 

2.81 

2.05 

2.66 

2.28 

1.92 

2.26 

2.09 

1.88 

4.0 

5.61 

5.30 

4.81 

5.24 

4.78 

4.66 

4.68 

4.62 

4.58 

4.5 

12.05 

11.97 

11.68 

11.74 

11.32 

11.57 

11.37 

11.16 

11.45 

5.0 

16.53 

16.09 

15.58 

16.19 

15.51 

15.28 

15.75 

15.38 

15.40 

6t (Gyr) 

1.0 

0.03 

0.02 

0.00 

0.02 

0.01 

0.00 

0.01 

0.01 

0.00 

2.0 

0.25 

0.18 

0.04 

0.14 

0.11 

0.01 

0.11 

0.07 

0.01 

3.0 

1.21 

0.92 

0.18 

0.79 

0.48 

0.02 

0.45 

0.31 

0.06 

3.5 

1.53 

1.17 

0.23 

1.06 

0.62 

0.02 

0.63 

0.43 

0.07 

4.0 

1.68 

1.31 

0.28 

1.21 

0.71 

0.01 

0.75 

0.53 

0.09 

4.5 

1.58 

1.36 

0.29 

1.18 

0.70 

0.00 

0.78 

0.54 

0.07 



Fig. 7. Run of the total H content as a function of luminosity for 
the 0.53512,0.61940, and 0.73821 Mq cooling sequences (solid, 
dot-dashed and dotted lines, respectively) Also shown is the H 
content for the 0.53512 Mq white dwarf sequence when element 
difussion processes are neglected (dashed line). The metallicity 
of the progenitor is Z = 0.0001. 


sequence, except for the short-lived stages at very high lumi¬ 
nosities (not shown in the figure), H burning never represents 
an appreciable source of energy. This is quite different from the 
situation for the 0.53512 Mq sequence. Here, the H-burning en¬ 
ergy release resulting from the proton-proton chain becomes the 
main source of energy over large periods of time at advanced 
ages, providing until logf ~ 9.4 a contribution to the surface 
luminosity that is much larger than that given by the release of 
gravothermal energy. The impact of proton-proton burning on 
the resulting cooling curve (solid line) is apparent. Indeed, note 
that towards the end of proton-proton burning phase, there is a 
change of slope in the cooling curve, until crystallization sets in 


at the center of the white dwarf and the cooling process slows 
down (for surface luminosities smaller than log(L/LQ) ~ -4). 
This is a result of the release of latent heat and gravitational en¬ 
ergy caused by carbon-oxygen phase separation. In passing, we 
note that the contribution of the CNO bi-cycle to the surface lu¬ 
minosity is restricted to the first 10^ yr of the white dwarf evolu¬ 
tion. However, except at very high luminosities, its contribution 
is always much smaller than the release of gravothermal energy, 
which drives evolution during these first stages. 

As mentioned before, all the white dwarf sequences com¬ 
puted in this work have been obtained following the progenitor 
evolution throughout the entire thermally-pulsing AGB phase. 
That is, we have not forced our models to abandon early this 
phase in order to get white dwarf structures. Given the fact that 
the mass loss adopted is just a rough extrapolation of what is 
known at higher metallicities, the initial-to-final mass relation 
offers a good testbed for our sequences. Then, we believe it is 
worthwhile to show the initial-to-final mass relation resulting 
from our full evolutionary sequences. This is done in Fig. |9] 
where our resulting white dwarf masses as a funct ion of the ini¬ 
tial ma ss at ZAMS are compared with the results of lRenedo et al.l 
(1201 oh for progenitors with Z = 0.01, and Z = 0.001, and 
with the final mass obtained by IWeiss & Fergusoril (l2009h for 
Z = 0.0005. In addition, we include the mass of the H-free 
core at the first thermal pulse as a function of the initial mass 
of the pro genitor for the case Z - 0.0005. It should be men¬ 
tioned that IWeiss & FergusonI (l2009h consider overshooting dur¬ 
ing the thermally-pulsing AGB phase, which considerably re¬ 
duces the growth of the H-free core. Thus, in this case, the fi¬ 
nal mass is not expected to be very different from the mass 
of the H-free core at the first thermal pulse. With this consid¬ 
eration in mind it is worth noting the good agrement between 
our mass of the H-fr ee core at the first t herma l pulse with the 
final mass found by IWeiss & Fergusoril (l2009h . By comparing 
our two sequences with Z = 0.0005 it is clear that, in the ab¬ 
sence of overshooting in the thermally-pulsing AGB phase, as 
assumed in our study, the mass of the H-free core grows consid¬ 
erably during the thermally-pulsing AGB. Hence, in this case, 
the final mass will be strongly dependent on the assumed mass- 
loss rates. This explains the steep slope of the initial-final mass 
relationship we find for the sequences presented in this study, 
since smaller mass-loss rates are expected at lower metallicities. 
Nevertheless, it should be emphasized that theoretical initial-to- 
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Fig. 8. Time dependence of the various energy sources and sinks for our 0.53512 Mq and 0.82623 Mq white dwarf sequences 
resulting from progenitors with Z = 0.0001. It is depicted the surface luminosity, L* (solid line), the luminosities arising from the 
proton-proton chain Lpp (dashed line), and from the CNO bi-cycle Lcno (dotted line), the energy loss by neutrino emission Ly (dot 
dashed line) and the gravothermal (compression plus thermal) energy release Lg (long dotted line). The domain of the main energy 
source is also indicated. Time is measured in years since the moment when the star reaches the maximum effective temperature. 
The progenitors masses are 0.85 and 2.5 Mq respectively. 


hnal mass relationships can not be accurately predicted, as they 
depend on two badly known processes in the thermally-pulsing- 
AGB phase, i.e., winds and the occurrence of a third dredge up 
episode. 

3.1. The importance of mass loss and overshooting 

We have shown that the mass of H that is left on the white dwarf 
at the beginning of its cooling branch is independent of the oc¬ 
currence of mass loss during the planetary nebulae stage. Here, 
we explore the role played by mass loss during the thermally- 
pulsing AGB phase in hxing the hnal H content of the white 
dwarf, and hence the importance of residual H burning. We con¬ 
centrate on those white dwarf masses for which we hnd residual 
H burning to be the main source of white dwarf energy. To this 
end, we re-computed the full evolution of our 0.80, 0.90, 1.0, 
and 1.25 Mq progenitors with Z = 0.0001 from the hrst thermal 
pulse on the AGB all the way to the beginning of the white dwarf 
cooling track, by artihcially increasing the mass loss rates by a 
factor of two. This is done during the entire thermally-pulsing 
AGB, as well as during the subsequent post-AGB stage. The re¬ 
sulting mass of H at the beginning of the cooling branch pre¬ 
dicted by these numerical experiments is displayed in Fig. [10] 
with dashed line. These predictions are compared to the results 
obtained whith our standard Z = 0.0001 progenitors with no 
enhanced mass loss. As expected, increasing the mass-loss rate 
yields smaller white dwarf masses. But note that for a given 
white dwarf mass, the hnal mass of H is the same, irrespective of 
the adopted mass-loss rate. In fact, the agreement in the relation 
H mass versus white dwarf mass is surprisingly good. This im¬ 
portant experiment shows that the H content with which a white 
dwarf of a given stellar mass enters its cooling track is inde¬ 
pendent of the adopted mass loss rates during progenitor evo- 



Fig.9. Theoretical initial-to-hnal mass relations resulting 
from our sequence s. In addition, we display the results of 
iRenedo et al.l (1201 Oh for progenitors with Z = 0.01 and Z = 
0.001, the mass of the H free core at the hrst thermal pulse for 
our Z = 0.0005 sequences , and th e initial-to-hnal mass relation¬ 
ship of lWeiss & Fergus^ (l2009h for the same metallicity. 


lution. Hence, our predictions that stable nuclear burning is the 
main source of energy over most of evolution of low-mass white 
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^wd^^sun 

Fig. 10. Logarithm of the total mass of H (in solar units) at the 
point of maximum effective temperature of the cooling sequence 
versus the white dwarf mass. The dotted (solid) line shows the 
predictions from our Z = 0.0001 progenitors when overshooting 
during the thermally-pulsing AGB is considered (neglected). In 
addition, we show the results for the case of an enhanced mass 
loss rate during the thermally-pulsing AGB and the following 
early post-AGB phases (dashed line). 


dwarfs coming from low-metallicity progenitors are valid inde¬ 
pendently of the rate at which mass is lost during prior evolution. 

On the contrary, and in view of our previous discussion, 
the occurrence of third dredge-up during the thermally-pulsing 
AGB phase is expected to reduce the final H content, and thus 
the impact of residual nuclear burning. To be more quantitative, 
we have re-computed some of our sequences but now allowing 
overshooting from the beginning of the thermally-pulsing AGB 
phase. It is well known that at low metallicities overshooting fa¬ 
vors the occurrence of a third dredge-up epis ode also for very 
small stellar masses (IWeiss & Fergusonll2009l) . To this end, we 
follow the evolution of our 0.85, 0.90, 1.0 and 1.25 Mq progen¬ 
itor stars with Z = 0.0001 from the first thermal pulse includ¬ 
ing diffusive overshooting with the overshooting parameter set 
to f - 0.015 at each convective boundary - see iRenedo et alJ 
(l2010l) : IWeiss & FergusonI (l2009h for details. Now, all these se¬ 
quences experience appreciable third dredge-up, with the result 
that the C/O ratio is strongly increased (see Fig. |2]l. Note that the 
C/O ratio reaches values of about 3-6, in g ood agreement with 
the predictions of lWeiss & Fergus^ (l2009il for Z = 0.0005 pro¬ 
genitors. The carbon enrichment of the envelope by more than 
three orders of magnitude as a result of a third dredge-up episode 
eventually leads to a more intense CNO burning, which strongly 
reduces the final H content of the resulting white dwarf, as can 
be appreciated in Fig. [10] Note that a reduction by more than a 
factor of two in the mass of the final H layer is expected if third 
dredge-up takes place. This strongly reduces the impact of resid¬ 
ual H burning during the white dwarf regime. We conclude that 
the occurrence of third dredge-up during the thermally-pulsing 
AGB evolution of low-mass stars (due to for instance to the oc¬ 
currence of overshooting) prevents H burning from being a rel¬ 
evant energy source in cool white dwarfs resulting from low- 



Log (Teff) 

Fig. 11. Evolution of the residual H content (in solar units) in 
terms of the effective temperature for our 0.9 Mq progenitor 
with Z = 0.0001 during the thermally-pulsing and post-AGB 
phases. Our standard sequence is depicted with a dot-dashed 
line. The sequence including overshooting during the thermally- 
pulsing AGB phase is shown with a solid line. Two additional 
artificial sequences with no overshooting are included: one in 
which the carbon and oxygen enrichment expected from the 
third dredge-up is only considered in the envelope opacity com¬ 
putation (dotted line), and another one in which the abundances 
of carbon and oxygen have been altered according to the predic¬ 
tions of the third dredge-up, see text for details. 


metallicity progenitors. In this sense, since the third dredge-up 
is strongly reduced in the sequences we presented in this paper 
(due to our neglect of overshooting), the final H content of our 
white dwarf sequences and the consequent impact of residual H 
burning on white dwarf cooling times should be considered as 
upper limits. 

In order to obtain a better understanding of the impact of 
overshooting on the final H content, we have re-computed the 
evolution of the 0.9 Mq progenitor star with Z = 0.0001 from 
the first thermal pulse, this time neglecting overshooting (and 
hence with no occurrence of dredge-up) but artificially taking 
into account the carbon and oxygen enrichment expected from 
the third dredge-up only in the envelope opacity. The result is 
shown as a dotted line in Fig.[TT] in which the resulting evolution 
of the residual H content (in solar units) is depicted in terms of 
the effective temperature during the thermally-pulsing and post- 
AGB phases. Our standard sequence with no overshooting is de¬ 
picted using a dot-dashed line, together with the prediction of 
the same sequence including overshooting during the thermally- 
pulsing AGB phase (solid line). It is clear from this figure that 
most of the reduction of the H content that characterizes the se¬ 
quences with overshooting is not the result of an opacity effect. 
Another artificial sequence neglecting overshooting was com¬ 
puted for which the abundances of carbon and oxygen in the 
entire model envelope have been altered according to the predic¬ 
tions expected from the third dredge-up (dashed line). In doing 
this, we neglect any other consequence of overshooting on the 
evolution. Clearly, it is the chemical enrichment of the envelope 
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that is responsible for the thinner H envelopes, not as an induced 
opacity effect, but as an increase in the H burning via CNO. 

4. Summary and conclusions 

In view of the recent finding that residual H burning in low- 
mass white dwarfs resulting from Z = 0.0001 progenitors 
could be the main energy so urce over most of their evolution 
dMiller Bertolami et alj|201^ , in this paper we presented new 
white dwarf evolutionary sequences resulting from progenitor 
stars characterized by a wide range of low-metallicity values. 
These new sequences are thus appropriate for the study of low- 
metallicity stellar populations. We computed the full evolution 
of white dwarf stars by consistently evolving the progenitor stars 
through all the stellar phases, including the ZAMS, the red gi¬ 
ant branch, the core helium flash (when occurred), the stable 
core helium burning, the AGB phase and the entire thermally- 
pulsing and post-AGB phases. Four progenitor metallicities have 
been considered; Z = 0.001, 0.0005, 0.0001 and 0.00003. To 
our knowledge, this is the first set of fully evolutionary calcu¬ 
lations of white dwarfs resulting from low-metallicity progen¬ 
itors covering the relevant range of initial main sequence and, 
correspondingly, white dwarf masses. A relevant aspect of our 
sequences is that the masses of the H-rich envelopes and of the 
helium shells were obtained from evolutionary calculations, in¬ 
stead of using typical values and artihcial initial white dwarf 
models. In computing this new set of white dwarf sequences, 
overshooting was disregarded during the thermally-pulsing AGB 
phase. Consequently, third dredge-up is much less efficient than 
in sequences that include overshooting at all convective bound¬ 
aries. Independently of the fact that evidence is not completely 
clear about the possible occurrence of appreciable overshooting 
during the thermally-pulsing AGB phase (even more in the low- 
metallicity regime), our reason for this choice is that we are in¬ 
terested in providing white dwarf sequences for which the resid¬ 
ual H burning has its maximum impact on white dwarf cooling 
times (see below). 

We found that in the absence of third dredge-up episodes, 
most of the evolution of low-mass white dwarfs resulting from 
low-metallicity progenitors in the range 0.00003 ^ Z 0.001 
is dominated by stable H burning. Specihcally, we hnd that for 
white dwarfs with M % 0.6 Mq energy released by H burning 
constitutes the main energy source at intermediate luminosities 
(-3.2 5 log(L/LQ) % -1), thus resulting in substantial delays 
in the cooling times. At intermediate luminosities, residual H 
burning increases the cooling times by more than a factor of 2 , 

and even at log(L/L 0 )-4, such delays amounts to 20 - 40% 

in low-mass white dwarfs. 

Our hnding that stable H burning is a main energy source 
that could be sustained over a long period of evolution in low- 
mass white dwarfs that result from low-metallicity progenitors 
opens the possibility that this stable burning could impact on the 
pulsational stability properties of these white dwarfs. The exis¬ 
tence of low-order g-mode destabilized by the e-mechanism in 
hot H-rich pre-white dwarf models has recently been reported 
dMaeda & Shibahashill2014ll . The assessment of fully nonadia- 
batic pulsation properties and a stability analysis of our models 
constitute an important issue, since this could provide for the 
hrst time evidence that e-destabilized modes could be expected 
in evolutiona ry models of cool, s tandar d carbon-oxygen white 
dwarfs - see ICorsico & Althaui (l2014t) for a recent report of 
evidence of e-destabilized modes in helium-core white dwarfs. 
However, this issue is beyond the scope of the present paper, and 
will be explored in a near future. 


We have also shown that our prediction that stable nuclear 
burning might be the main source of energy over most of evo¬ 
lution of low-mass white dwarfs resulting from low-metallicity 
progenitors remains valid independently of the rate at which 
stellar mass is lost during the thermally-pulsing and post-AGB 
phases, and during the planetary nebula stage. However, we find 
that the mass of H of the resulting white dwarfs becomes sustan- 
tially smaller (by a factor of two) as a result of the occurrence 
of third dredge-up during the thermally-pulsing AGB phase. In 
this sense, we hnd that the carbon enrichment of the progenitor 
envelope resulting from the third dredge-up prevents H burning 
from being a relevant energy source in cool white dwarfs result¬ 
ing from low-metallicity progenitors. 

According to our simulations, white dwarfs populating low- 
metallicity globular clusters could be characterized by signif¬ 
icant residual H burning, so as to leave some signature on 
the observed white dwarf luminosity function of such clusters. 
NGC 6397, a metal poor ([Fe/H]= -2.1), old globular clus¬ 
ter, whose well-determined white dwarf luminosity function 
has been proposed as a tool to constrain white dwarf physics 
(IWinget et al.ll200^ is particularly relevant in this context. In 
particular, this cluster could be used as a testbed for the occur¬ 
rence of stable nuclear burning in their white dwarfs, which, in 
view of the results of our study, could be used in turn to con¬ 
strain the efficiency of third dredge-up episodes and the ensuing 
carbon enrichment during the thermally-pulsing AGB phase of 
low-metallicity progenitors. 
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